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Abstract 

We consider the quantum mechanics of an harmonic oscillator when it is driven either by an 
external random (white noise) force or when its frequency is sinusoidally time-dependent, either 
varying slowly (adiabatically) or at about twice the natural frequency (parametrically). We focus 
on finding transition probabilities and expectations, for which it proves convenient to utilize the 
Wigner-Weyl formulation associating an operator A with its phase space transform, A •<—> A(p, q) 
and in particular the density matrix p with its corresponding Wigner function defined by p(t) < — > 
hpwiPjq'A). It is shown that for long times, whatever the oscillator’s initial state, a white noise 
driving force randomizes over time the distribution with respect to phase for any operator $ —> 
$(<()) and in particular for the phase operator cj) <—> cj) itself, where cj) is the polar angle in the phase 
plane. We also consider, in the long-time limit, the randomizing effect of the white noise force for 
expectations of (jP and of any operator Cl < — > B(i?), where R is the radial coordinate in the phase 
plane, the partner to (p such that ^ -I- When there is no external force acting, 

but there is a time-dependent frequency w^(t) = ujq{ 1 + e(t)), where e is small, then, if e varies 
slowly, transition probabilities are unaffected, but for a parametric oscillator e varies sinusoidally at 
a frequency of about twice uig and these probabilities are strongly time-dependent. 

PACS: 03.65.-w; 03.65.Ge; 05.40.-a; 42.50.Lc 


1 Introduction 

Analyzing the quantum harmonic oscillator—or more generally any polynomial Hamiltonian quadratic 
in p and q with time-dependent coefficients—is about as old as quantum mechanics itself. A recent 
approach [T] summarizes and develops an analysis based on dynamical invariants [2] to calculate 
the time evolution of wave functions. One can also approach the problem using path summation [3]. A 
third alternative (at least when the frequency and/or mass depend on time) is solution by ansatz [4]. 
These methods can be, and have been, used to calculate the effects on its wave function of a given time- 
dependence in an oscillator’s Hamiltonian. In the present paper, however, the problem is transformed 
to phase space using the Wigner-Weyl association. It appears that this method is particularly suited 
to calculating time-dependent expectations and transition probabilities, and that is our focus here. 

There are of course models of the quantum mechanical response of a small system such as an 
oscillator to the action of a heat bath [5] or also with with external forces acting [6] . Here we analyze 
exactly the response of certain operators associated with an oscillator, to the randomizing effects of 
a white noise random external force, for which the transformation to the Wigner-Weyl picture proves 
useful. 
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In general, the association of an operator, A <—?• A(p, q) with a function on the phase plane can 
take many forms [7], but the commonest, and in many ways the simplest, is that of Wigner and of Weyl 
HlElllo]. This paper utilizes only the Wigner-Weyl correspondence. The density matrix p and phase 
operator cj) (see reference m) are well-defined quantum mechanical objects with the correspondences 
Pw{p,Q',t) <— > p{l)/h and — >(f>, where pw{p,q]t) is the time-dependent Wigner function and (p 

is effectively the polar angle in the phase plane {p, q). Here I mostly consider the time dependence of 
operators related to the one-dimensional harmonic oscillator under the influence of a time-dependent 
force or when its frequency is time-dependent. The general Hamiltonian for these cases is 

- Kt)q ■ ( 1 ) 

A great advantage of using the Wigner-Weyl phase space formulation for the oscillator is that when 
the quantum mechanics is ‘unwrapped’ the dynamics is purely classical. A price to pay, however, is 
that to make the final result properly quantum mechanical one must perform one or more integrals 
over phase space. 

Section [2] reviews briefly the basic association A i—> A{p,q). The time-dependence of A and its 
application to the harmonic oscillator, equation ([I]), is discussed section [3l The material in sections 
[2] and [3] is not new. In section H] we consider the action of a white noise Gaussian random force A(t) 
whose ensemble average vanishes but with X(ti)X{t 2 ) = p6{ti — 12 ). We use an analysis of classical 
Brownian motion in phase space |12] to find the ensemble averaged quantum mechanical propagator 
for the Wigner function. From this one can calculate transition probabilities for the oscillator. In 
particular we calculate explicitly the probability that the oscillator leaves the ground state after a time 
t under the action of X{t). 

Section [5] is devoted in the first instance to the time evolution of the phase angle operator cp 
and then, more generally, to operators $ < — ^{p) under the influence of the stochastic force. An 
expression is derived for the expectation Tr(/5(t)$) when the initial state is the oscillator’s ground 
state |/io)(ho|. More generally it is also shown that, whatever the initial state may be, at long times, 
the distribution of ^{p) is randomized with respect to the angle p. Section [5] also considers the long 
time effects (for any initial state /o(0)) of the random force on the expectation Tr(p{t)Cl) of operators 
H —)• Xl{R) where R is defined by the Weyl transform of the free oscillator’s Hamiltonian such that 
p^/{2m) + /2 = hjjR?/2. Section [5] ends by considering, in the long time limit, the influence of 

the random force on the behaviour of the expectation TXc^piipp'^), for any initial state p(0). Although 
p"^ itself has a Weyl transform depending on both p and R, this expectation becomes randomized in 
time with respect to p. 

Section [6] considers, approximately, two special cases of a time-dependent frequency uj{t), but with 
no other external influence. On the one hand when is slowly varying (the adiabatic limit) a 
simple time dependent phase shift results, generally of limited physical consequence. However, for a 
parametric oscillator, the frequency varies sinusoidally at about twice the oscillator’s natural frequency. 
In that case it is shown that when the oscillator is initially in the ground state the probability of it 
staying there decreases, in the limit, exponentially in time. 

Section [7] gives a brief discussion. 


2 The Wigner-Weyl picture 

There are many possible formulations of quantum mechanics in phase space [7]. Generally, they can 
be related [Hill] to that of Wigner and Weyl, for which we shall adopt a formal but efficient notation 
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m- Denoting the Weyl transform of an operator A by A{p,q), or sometimes by {A)(p,q), it is given 
by 

A <—^ Tr(i A{p, q)) = A{p, q) = {A){p, q) 


where 


A{p, q) = 




. ^ z 


and D is the Weyl operator l[T6]b 

b{p,q) = 

Formally A has the the properties m that its trace is unity, that 

dpdg 


f 


h 


-A{p,q) = 1, 


and that 


Tr 


(a(p, q)A{p', g')) = h5{p - p')5{q - q’) . 


Using this, equation (l2|) can be inverted to give 


A = 


dpdq 


h 


A{p,q)A{p,q) and TtA = 


dpdq 

h 


A{p, q) ■ 


From these properties one can also show that, for two operators A and B, 

' dpdq 


Tv{AB) = 


h 


-A{p,q)B{p,q) 


and that the Weyl transform of the product AB is 


AB <—Tr 


(^ABA{p,q)'j = A{p,q) 


exp 


\h d d* d\ 

2 \dq dp dp dq) 


B{p,q) 


( 2 ) 


(3) 


(4) 

(5) 

( 6 ) 

(7) 


( 8 ) 


(9) 


where the starred operators act to the left on A{p, q). 

The Wigner-Weyl association can be expressed in terms of either D or A, whichever is more 
convenient, for they are related by the first of equations ([3]), or equivalently, by 


D{p, <l) = j 


dpdq , , 


■ e ft 


A(p',g'). 


( 10 ) 


Algebraically, the relation between D and A can be expressed [Ullla] 


as 


A(p, q) = 2 D{2p, 2q) H, so that D{p, q) = -^ A(p/2, q/2) H, 

where 11 is the parity operator. 

The Weyl transform makes the fundamental association 

b{p',q') = eH(p'5-9'P), 


( 11 ) 
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so that, sensibly, — >1, F{p)i — >F{p) and F{q)i — >F{q). But functions that mix p and q are more 
complicated. For instance p q< — >pq—\h/2 and q^pi — >p q+\h/2 so that the Weyl operator corresponding 
to pg' is {pq + qp)/2. The Weyl transform of the Hamiltonian operator, equation ([1]), takes the classical 
form ^ 

H{t) i —^ + Kn(^‘^{t)q^ - X{t)q (12) 

Finally, the the Wigner function is defined as 

Pw{p, q; t) = ■^Tr(p(t)A(p, q)) < —> ^p{t) , (13) 

so that by ([8]), 

/ CXD 

dpdqp^,{p,q]t)A{p,q). (14) 

-CO 


3 Time development 

The evolution of the Wigner function under the action of a time-dependent Hamiltonian is is well-known 
[laiaiis]. Wave functions evolve according to 

\i^t) = Utl'p), (15) 

where the unitary time evolution operator Ut is governed by the equation 

ih^Ut = HtUt, (16) 

and we allow for explicit time-dependence in the Hamiltonian Ht. Then the time-dependent density 
matrix p{t) is given by 

m = utpmui, (17) 

and its Weyl transform is 

Pw{p,q-,t) = j dp'dq'P^{p,q,t\p',q' ,0)Pw{p' ,q'-,0) (18) 

where Pw{\) is the Wigner propagator dehned by 

Pw{p,q,t\p',q',0) = ^Tr{U^A{p,q)UtA{p',q')). (19) 

In particular, it is easy to show from the properties in section [2] that 

J dp'dq'P^{p,q,t\p,q',0) = J dpdqP^,{p,q,t\p',q',0) = 1 (20) 

and 

Pw{p, q, 0|p', q', 0) = 6{p - p')6{q - q'). (21) 

From dehnition (fT^ and equation (fTHl) we can differentiate this propagator with respect to time to get 
its equation of motion as the Weyl transform of a commutator, as follows: 

O -I 

^f^-^Pwip',q',t\p,q,o) = j^[[Ht,UtMp,q)Ut]] ip'^q') ( 22 ) 


4 









The Weyl transform of Ht given by equation (jl2p is quadratic in p and q and so has no derivatives with 
respect to p and/or q of order higher than second. Applying equation ([9]) to (l2^ is straightforward: 
collecting terms gives 


■^Pwip,q,t\p',q',0) = 


(23) 


This describes classical motion under the action of Ht{p,q). The solution that satishes condition pill 
is 

Pw{p,q,t\p,q,0) = ^{p-p{t\p',q',0))S{q - q{t\p,q',0)), (24) 

where {p{t\p', q',0), q{t\p', q', 0)) is the classical phase space solution for momentum and position under 
the action of Hamiltonian Ht{p,q) such that {p{0\p', q', 0), q{0\p', q', 0)) = {p',q')- This solution also 
obeys equation (|20p . by direct integration in the hrst instance and, in the second, by recognizing that 
the Jacobian d{p, q)/d{p', q') is unity when {p, q) and {p', q') are related by the classical motion implied 
by equation (|24ll . For the cases we are considering, the equations for classical motion are 


p{t) = mq{t) and p{t) = —muP‘q{t) — A . 


(25) 


where w and/or A may depend on time. 


4 Forcing by stationary white noise 

4.1 Time dependent external force 

As shown in section [3l the time evolution for the harmonic oscillator can be found by solving the 
classical equations of motion. For instance, when the oscillator is subjected to a time-dependent force 
only, then A depends on time and oj and m are constants. Expressed in terms of dimensionless variables 
(x,y) = {p/{ha), aq), where = mui/h, equations (1^ can be written 


^(x + iy) - iuj{x + iy) = 


(26) 


Writing x{t) and y{t) for x(t|xo, yoi 0) and y{t\xo,yo,0), the solutions with initial values xq and yo, we 
have 

x{t)+iy{t) = {xo + iyo)-\[ dse"‘‘^^A(s) . (27) 


ah , 


For example, when 


where 


X{t) = Ao sm{ujt + 0 ), 

x{t) + iy(t) = - sin(a;t)e-‘^) , 

2ahu} V / 

xo + m = RoP^°- 
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4.2 The motion 


We model the effects of a stationary random force by supposing the Hamiltonian has the form 




+ -n-Q - Ht)Q, 


(28) 


2m 2 

where X{t) is a white noise stationary Gaussian process [19] with the particular ensemble averages, 

A(t) = 0 and X{ti)X{t 2 ) = ^S{ti — t 2 ) ■ (29) 


One can, of course, use equation (l271l and take the relevant ensemble averages, but this has already 
been done for the classical theory of Brownian motion m- In that theory, a traditional approach is 
through the Langevin equations. In one dimension these are 


p = mq and p = F{p, q) — Pp + X{t ), (30) 

where X{t) is the random force and F may depend on p or q. The friction term, —/3p represents the 
effect of a heat bath and cannot be modelled by any Hermitian single particle Hamiltonian. Here we 
shall be considering the limit /3 —)• 0. Denoting the the ensemble average of the classical propagator 
for Brownian motion by an overline, from equation (1241) we have 

Pw{p,q,t\po,qo,0) = S{p-p{t\po,qo,0))S{q - q{t\po,qo,0)) = W{p,q,t\po,qo,0) 


where W{-\-) is the conditional probability density for [p, q) at time t given the initial conditions {po, qo) 
at t = 0. 

Classical Brownian motion is random, stationary and Markovian, and it is characterized by the 
conditional probability density W{p, q, t\po, qo,0). Reference [T2| gives the following efficient expression 
that is generally approximate, but exact when F{p,q) is at most linear in p and/or q: 


W{p,q,t\po,qo,0) ~ 


dad6 


exp(io(g - g(t))exp(i6(p - p(t)) 


X exp 


“2 / + 


(31) 


where p characterizes the random force X{t), {pit),q{t)) is shorthand for the solutions 
(^(^Ipoj ? 0 ) 0), g(t|p 0 ) 90 ) 0)) to (l30|) (most generally with /3 / 0), and dp^ is the partial derivative with 
respect to the momentum ps at the intermediate time 0 < s < t. Thus p occurs explicitly in equation 
()31l) and /? occurs implicitly through the solution {p{t), q{t)). If the classical Brownian particle were to 
thermalize after long times |I2] then the relation ^/2 = l3kT/m must obtain, where k is Boltzmann’s 
constant and T the temperature. Here we choose to consider the limit /3 —>■ 0, so modelling an oscillator 
driven only by the random force A(t), equations (I28h and (I29p . Another way to look at this model (see 
reference |6]) is via the WFPE (Wigner-Fokker-Planck equation) for an oscillator in thermal equilibrium 
with a thermal heat bath in the limit of high temperatures but small friction such that the product 
f5kT remains finite. 

In our picture we require a solution for the unforced harmonic oscillator, p{t) and q{t), with van¬ 
ishing (5. These are, for 0 < s < t, 


and 


p{t) = Ps cos {ui{t — 
q{t) = 

moj 


s)) — muqs sin (a;(t — s)) 

(32) 

- s)) + qs cos{u}{t - s)). 

(33) 
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Using this information to evaluate (j3ip and replacing ~ by =, for this result is exact, gives 

dad6 


W{p,q,t\pQ,qQ,Q) = 


. (2^)2 

X exp(i6(p — p(t))exp 


exp(ia(g - q{t)) 


— — / ds sin (US + 6 costas') 
2 Jo Vm(u / 


(34) 


wherein p{t) and q{t) refer to the solutions with initial conditions pq, q^. The simplest case is the free 
particle, for which (u ^ 0. Taking that limit in (I34p gives, in detail. 


Wip,q,t\po,qo,0) = 
f dadb 


(2vr) 


exp (^ia(g - go - exp(i6(p-po)) exp (~'® + ^) 


(35) 


All the integrals in and (f35p can be easily evaluated if necessary. 


4.3 Transition probabilities 

The probability for transition between states {tpi) and \^p 2 ) is 

= j J dp'dq'{\i’ 2 ){ 1 ^ 21 ) {p,q)Pwip,q, t\p',q',0) {p',q'), (36) 

where is the Wigner propagator, (jl9|) . In particular, for the random driving force p29|) . the ensemble 
averaged propagator is 

W(p, q, t\p', q', 0) = P^{p, q, t\p', g', 0). 

The corresponding ensemble averaged transition probability is 

|(V'2|Ut|V’i)P = j j dp'dg'(|V'2)(V’2|)(p,g)lU(p,g,t|p',g',0)(|V’i)(V'i|)(p',9') (37) 

where, for an oscillator, W is given by p4P . Writing {dx,dy) for partial derivatives with respect to x 
and y, where {x,y) = {p/{ha),aq), and o? = muj/h, that expression for W can be formally rewritten 
as 


W{p,q,t\p ,t) = exp 


N r*' 

— / d0 (sin0cly + cosScIa;)^ 
2 Jo 


6{p - p{t))5{q - q{t)). 


with 


6{p — p{t))6{q — q{t)) = ^ S{x — x{t))S{y — y{t)) and the definition N = — 

n muj^n 


(38) 


(39) 


In (|39l) {p{t),q{t)) is the solution for the free oscillator with initial conditions {p',q'), that is to say 
p{t) = p{t\p', q' , 0) and q{t) = q{t\p', q', 0), equations (1521) and (I55P . so this classical propagator conserves 
the oscillator’s energy. 

As an example let us suppose that the oscillator is initially in the ground state, \ho), and seek the 
probability in time that it stays there so that IV’i) = IV’ 2 ) = l^o); where 


(Cl^o) = (Q!/\/7r)^^^exp (—q;^^^/ 2) and . 

The Weyl transform of the initial (and final) state is easily found to be 

{\ho){ho\) {p,q) = 2exp {-{x^ +y^)) = 2ex.p{-R^). 


(40) 


(41) 


7 















This sole dependence on R suggests that it will be convenient sometimes to integrate with respect to 
plane polar coordinates (-R, (j)). This can be effected in the integrals by the equivalences 

dpdq = hdx dy = h R dR dip , 


and the free motion in these coordinates is 


6{p - p{t))6{q - q{t)) = - R')5{(t) -cp' - ujt) 


Therefore, performing the integrals over R! and (p', gives 

dpdq 


\{ho\Ut\hoW = J 

= 4 


h J 

dpdq 

h 


dp'dq' {\ho){ho\) {p, q)W{p, q, t\p', q', 0) {\ho){ho\) {p', q') 


exp(—exp 


N 

— / d0 (sin05y + cos^^a;)^ 
2 Jo 


exp(— 


We can now make the Fourier decomposition 


exp(-i?2) = J +^2 ^exp(i(fcix + fc2i/)) 


so leading to the expression 
dx dy 


\{ho\Ut\h^)\^ = d 


2tt 


exp(— 


dkidk2 


dvr 


exp 


t-2 ^2\ 

\ ^ ]exp(i(kix + k 2 y)) 


X exp 


jY rujt 

'Y 


riJt 

/ d9 {ki cos 6 + k 2 sin 9)^ 
Jo 


(42) 


(43) 


(44) 


Expression (1441) is a four-fold quadratic integral and thus can be evaluated exactly. We can start with 
the double integral with respect to ki and k 2 - The standard Gaussian form for this is 



dx"'exp 



Ax 


exp[i 


X 


(27r)’^ 

detA 


exp 



(45) 


where A is a real symmetric n by n matrix. This sort of integral finds application, for example, in 
evaluating path sums in field theory [20], often in the limit of infinite n In the present case, with n = 2, 
= (/ci, A: 2 ), and = (x, y), we have 



f5l2 ^ 

where S = 

^ d9 cos^ 9 d9 sin 29 ^ 

1 f^l2 

^+A522 j 


d9 sin 29 d9 sin^ 9 j 


The inverse of A is 


A-i 


1 

detA 


Combining (1440 and (1451) then gives 


\+NS22 -f5i2 \ 

-f5l2 \+NSll y 


\{homhoW 


ttV detA 


dx dy exp (—) exp 



(46) 


(47) 


(48) 






























This has the standard form 


dx ^exp 


-X 

2 


B X 


(27r)»^ 

detB 


where, in this case, n = 2 and x = (x, y). For ()19]) it is easy to see that 


B = A’l + 21 


(49) 


where I is the unit matrix. Finally then, 

_ 2 

\{ho\Ut\ho)\ - ^et(2I + A-i) “ det(I + 2A) 

= , ^ , (50) 

^(l + ^f-(f)%in2.t 

where N is given by (1391) and we have expanded the determinant using expressions ()46F 

5 Phase 

5.1 Definition 

Generally in the Wigner-Weyl picture, the time-dependent average of an operator A with respect to 
the state p is given by (fT4|) with (fT8|) . where, for an oscillator, the evolution is classical, equation (1241) . 
In particular consider the the Weyl quantization (jx — >(j)(p,q), of the harmonic oscillator phase in the 
plane (p, q). Properties of cj) have been considered previously [HI |2T1 |22l |23] . It is a bone-fide bounded 
self-adjoint operator on Hilbert Space. As befits an angle, its spectrum must be limited to a range of 
27r, but it appears that a published proof mi [ 22 ] only exists limiting the spectrum to a range of Stt. 
There may, however, be an unpublished proof [24|— using the ‘method of wedges’ defined in section 
9.4.5 of [H] — limiting the spectrum to a range of 27r. In any case, numerical calculations representing (p 
by the finite matrix {hm\(p\hn), 0 < {m,n) < N in terms of the harmonic oscillator energy eigenstates, 
indicate—as N is increased up to several hundred—a smooth convergence from within to a uniform 
spread of eigenvalues over an interval of length 27r. In the following we consider the time-dependence 
of 4> when the oscillator is driven by a white noise stochastic force, equations (1281) and (1291) . 

We choose to define 4> in terms of the harmonic oscillator creation operator 

at ^ ^ (“« - «) = (“« - ‘is) 

where, for an oscillator having mass m, = muj/h . Then, in terms of the dimensionless variables 
(x, y) = (^, aq) the Weyl transform of af is where E? = x'^ + and the the energy of a free 

oscillator is ^ ^ The Weyl transform of (p involves an integral over 

the phase plane {p,q). So if we use plane polar coordinates {E,(p), where dpdq/h = i?di?d</>/27r, we 
must restrict (p to a range of 27r which we shall take to be the interval [—tt, vr). Then 

(p = TT{$A{E,(p))i—= dEE j ^(pA{E,(p), (51) 

where we now consider A as a function of plane polar coordinates E and (p so that the matrix elements 
of (p with respect to the energy eigen states |/i„), (n > 0) of the harmonic oscillator, with energies 
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(n + l/2)huj, are given by 


where [22] 


{hm\^\hn) = j dRuJ ^ (j) {hm\ (p) \hn) , 


(52) 





(53) 


Here n(,{ng) is the lessor (greater) of the pair {m,n), and is the Laguerre polynomial [25]. When 
(l53]) is used in ([5^ one finds after some manipulation [HI [22] that 

•n—m+1 

{hm\ (plhn) — (1 Sm,n) 9m,n 

m — n 

where gm,n is the symmetric matrix 


_ r (^ + s^) 

~ r (!^ + s,) V nel 


(55) 


with 


se = 


1 /2 n£ even 
1 ni odd 


(56) 


Alternatives to (j) have been suggested to represent the phase of an harmonic oscillator, or in 
some sense photons, for instance by a non-projective positive operator valued measure, or POVM 
[MiisaiMiissiisniiM]- Using a POVM to represent a quantum system may be thought of as allowing 
for an element of imperfection in the measurement. 


5.2 Forcing of angle functions by white noise 

Here, let us consider those operators, $, whose Weyl quantizations are functions only of angle, 
in the phase plane, thus 

—t<h(())) 

Then, let us ask for the expectation of such a $ when the oscillator is initially in the ground state, ho, 
and subsequently driven by stochastic Hamiltonian (|28p with ()29l) . namely 

Ti{p{t)^) = J J dp'dq' W{p,q,t\p',q',0) p^{p',q',0), (57) 


where 

Pw{p,q,0) = {\ho){ho\) {p,q) 

is the Weyl transform of the initial state. The analysis of equation (I57p is similar to that used on 
equation ([37]). The result is 


Tv{pm) = 


1 r d(j) 

VdetA J-TT 2vr 

r 


roo 

^{(p) / di?i?exp 

Jo 


27r 


<^>{cP) 


j^2 , . 

——— ( 1 + Vwt — N sin Lot cos(2o!) — Lot ')) 
4detA V \ -r jj 

a/(1 + Niot)'^ - A2 sin^ cot r dcp 

= J —^{P))P{<P,ujt), 


[(1 + Not) — N sin iot cos{2(p — ot)] 


(58) 
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where we have expanded 


(59) 


4 det A = (1 + Nuit)"^ — iV^ sin^ ojt. 

Figure 1 shows a computed plot for the choice 4’(0) = 4>. 

The factor P{(j),ujt)/{27r) in (IHHI) is a bone fide probability density. To see this, first note that, for 
all t > 0, it is positive. Secondly, when sinwt vanishes, and also in the limit t —)• oo, it reverts to l/27r, 
and so corresponds to a random distribution of phase and is manifestly normalized with respect to the 
interval [—7r,7r). When sin cut does not vanish we can, after some algebra, re-express the integral (I58p 
(with 4>(^) = 1) in terms of the variable z = exp(2i(p) as 





—sign(sincut) exp(icut) / 

Jc 


dz — 1 

27ri {z — z^){z — z~) 


(60) 


where the curve C is taken twice anti-clockwise around the unit circle in the z-plane, and 


z^ = exp(icut)(A ± VA^ — 1) with A 


1 -|- Noot 
N sin Ldt 


Finally, the integral in (I60p can be seen to be unity, for when sign(sinc<;t) = -|-1, z"*" lies within the unit 
circle and z“ lies without, but when sign(sincut) = — 1 these roles are exchanged. 


5.3 Forcing by noise at long times 


Equation (f58ll . an exact result when the initial state is \ho), indicates that the random process X{t) 
randomizes the phase phase distribution after long times. One might reasonably expect this to occur 
whatever the initial state, pw{p,q,0). To verify this, from equation (1571) it is clear that we need the 
limiting behavior of W{p, q, t\p', q\ 0) for large t. 

The general expression for VF(|) is given by (j34l) . For long times we can make the approximation 


[ ds f-sincus-|-^coscus) ~ 

Jo / 


(-) ^'-^ 


In this limit, the integrals over a and b in (|34l) reduce to the product of two simple Gaussian integrals. 
Expressed in terms of the dimensionless variables x and y the long-time limiting behavior of W is 


W{x,y,t\x',y',0) ~ 




1 

Nujt 


{x 


Xtf + {y- ytf , 


where xt and yt are the oscillator’s free motion with initial values x' and y' at time zero. In polar 
coordinates we can write 


X + \ y = R expi(/) and xt + iyt = R' expi(())' + ojt ), 

so that 

[x — xtf' + [y — yt)^ = R? + R'^ — 2RR! cos{(j) — JJ — ojt). 
Then using this in (I57p gives (for long times) for any operator 

Oi — >0{p, q ), 


Tr(p(t)0) 

1 

^ Niot Jq 



dRR 0{R,4>) 


1 

Nujt 


{R^ + R’^ 


2RR! cos((/> — — ujt)) 




(61) 
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In particular, for those operators <I> whose Weyl equivalents <&(</>) are functions only of (j) we can write 
(changing variables to a: = R/y/Nojt) 


Tr(p(t)<h) 


dcf) 


poo p 

I 


— TV 

TT 


m f 

Jo 


dx X 


-exp 


, 2 R xR' , ^ 

— 1 X + —7-C0S((^ — (j) — Ojt) 


Nut y/Nut 


PniiR'A^,0) 


dcj) 


/■oo roo pn jil 

<!>((/>) / dxx / dR'R' / —^exp (—x^) (/>', 0) 

Jo Jo J-TT ^ 


dR'R' ^p^(R'^cj)',0)= ^$(</,)Tr(p(0)) 


Thus for these operators the phase distribution is randomized in time by the stochastic force. 

Similarly, we may ask for the long-time effect on operators whose Weyl equivalents are functions 
only of R, namely those n(i?) such that 

n< — >n{R). 

In this case, the same sort of formal manipulations give 


Tr(/5(t)ll) — y [ — [ dx X Q{xV Nut) 

J-^ 27t Jq 

.TiR'R' r Me: 

Jo J—n ^ 


-[x^ + 


xR' 


: COs{(j) — (j)' — Ut) 


Nut y/Nut 

/ TT fOO 

— / dxx n(x\/ Nut) 2 exp(—x^) Tr(p(0)) 

-TT 2vr Jq 
poo 

/ dx exp(—x) n(\/X Nut). 

Jo 


Pu,iR',(j)',0) 


(62) 


For instance if we choose 11 = exp(—where /3 is a parameter and H = (h + ll2)fiu is the 
Hamiltonian for a free oscillator, then, from equation ()53jl and reference to standard tables (eg [25]), 


exp(—/3i7) Tr ^exp(—/3iJ)A(i2 ,1 

CX> 

= 2 ^ exp( — {m + 1/2)1tu(3) {—)'^exp{—R^)Lm{2R^) 


m=0 


cosh ( ^ 


exp 


—R^ tanh 


hul3\ 

^ J 


Thus, in this instance we have, using (I62p and considering the limiting behavior for large t, 


Tr 


i- (^p{t)exp{-/3Hfj 


cosh 


/ hufJ 


+ Nut sinh 


/ hufJ 


-1 


(63) 


Although the Weyl equivalent to (j) is the pure phase 4>, the Weyl equivalents of powers and functions 
of cj) are often not functions of (j) only. Indeed, even the Weyl equivalent of the square (fP' is a mixed 
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(64) 


function of 4> and R. One can show ([ 21 ]) using (j54p and choosing the range [—7r,7r), that 


LXJ IIL LX; 

{hm\i’ \hm) = ^ ^ ^ ^ {9m,m—n) + ^ ^ {9m+n,m) j 

n=0 n=l n=l 

where the first sum on the right-hand side only contributes when re > 1. It is, however, a property 
of 9m, n, equation (|55]l . that as both m and re increase to values that are large with respect to their 
difference (the correspondence limit) then gm,n 1- Thus as m increases in (I64h the standard deviation 
{hmw\hm) approaches 7r^/3, a value associated with a random distribution of phase ([2T]). 

We may ask what effect the random force \{t) might have on In particular, if its Weyl transform, 
in polar coordinates, is 

= , (65) 


( 66 ) 


then, adapting the development leading to (f62]l . we have (for long times) 

>2 f dxxexp(—x^) [ — (cp) (xV Ncot, 4>). 
Jo J- 7 T 27r 


TripitW) 

Now, from ([53]) and (l65]) . 


/ •yj- 1 I tXJ y 2 2 \ 

■£{^'^){RA) = + y 2(-)™exp(-i?2)L^(2i?2) 

-TT ^ / 

OO y 2 \ 

+ ( {hm\^^\hm) - ^ ) 2{-reM-R^)Lmi2R ^), 

771=0 ^ ^ 


771=0 

9 OO 

TT 

T 


where I have used (in the limit z ^ —1) the equation 


(1 — z) ^exp 


xz 

z-1 


Lrn{x)z^‘ 


771=0 


Using this in equation (jM|) shows that for long times 


Tr(p(f)(/>2) —^ — . 

This is characteristic of a random distribution of (j) over the interval [—vr, vr), namely 


/: 


2n 


TT 

Y 


(67) 


6 Time dependent frequency 

6.1 Adiabatic frequency dependence 

We apply the ‘method of averaging’ [32| to the oscillator, equations (12511 when there is no external 
force but with a time-dependent frequency, 

Y(f) = a;g(l-be(t)) , (68) 

where wq is a constant base frequency and we will ultimately assume that e(f) is very small compared 
to unity. Then defining variables (x, y) = {p/(ha), aq), where = mojo/h, we may opt to express the 
(unknown) solution to the equations of motion in the form 

x(t)+iy(t) =r(t)e'(‘‘^°*+®W). (69) 
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Then 


x{t) + iy{t)= r{t)+ ir{t){uJo + 9{t)) . 

This equation is a result of the choice (169 j) . But the motion must be governed by the equations 

p{t) = mq{t) and p{t) = . 


(70) 


These require that 

x{t) +iy{t) = ^^ + iaq{t) = -^u^{t)q{t) + i^p{t ), 
which can rearranged to give 


x{t) + iy{t) = iujQ{x{t)+ \y{t)+ ie{t)y{t)) 

= ia;or(t)ei(-«*+^W)-woe(i)y(t) (71) 

Equating the right-hand sides of (l7Up with (I69|l . and (EH), simplifying, and then separating the real 
and imaginary parts gives 


r{t) = -^ujor{t)e{t)sm{2{uJot + 9{t)) 
9{t) = Wo e(t) sin^(wot-b 0(t)). 


(72) 


For the adiabatic limit we assume that base period, is very much smaller than the time over 

which e{t) varies signihcantly. In that case, denoting by 9a{t) and ra{t) a ‘local’ or adiabatic time 
average of 9 and r over a time long compared to the base period 27r/a;o but short compared to e(t)/e(t) 
we have ^ 

ra{t) = 0 and 9a{t) = - woe(t), 

so that 


^a(t) = ^(0) and 

(73) 

9a{t) = 9{0) + e{t), where 

(74) 

e{t) = Y dse(s). 

(75) 


To this level of approximation the simple time dependent phase shift e{t) will not affect transition 
probabilities. 


6.2 Parametric oscillator 

Suppose the frequency varies at approximately twice the fundamental rate wq, so that 

a;^(t) = Wq)!- becos (2(wo-b/)t)], (76) 

and that e and //wq are small compared to unity. Then, assuming an approximate solution in the 
form of equation (|69]) , and defining the parameter u = ewo /4, equations (|7^ take the form 

f{t) = —2ur{t)cos{{2u}o + f)t)sin(2{LOot — ft)) 

9{t) = Aucos{{2loo + f)t)sm^{u}ot + 9{t)). 
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Expanding the right-hand sides of these equations using basic trigonometric identities, and retaining 
only the slowly varying terms gives 


r{t) ~ —ur{t)sm{29{t)—ft) 

9{t) ~ —u cos{29{t) — ft). 

This model is discussed in reference [33], where it is shown that growth occurs in the classical motion 
for values of / such that —2u < f < 2u. For simplicity we choose to pursue the case / = 0. Then the 
equations to solve are 


r ~ —ur{t)sm{29{t)) 
9{t) ~ —ucos{29{t)). 

The second of these can be solved by direct integration: 

tan (^9{t) + = exp(—2ut) tan 


(77) 

(78) 


(79) 


where 0(0) = (j)Q is the initial angle in the phase plane for the oscillator’s motion. Consistent with this 
we assume the following solution for x{t) + iy{t): 

(80) 


x{t) + iy{t) = Rq exp (^i(ti;ot “ ^) j exp(ut) cos ^ j + iexp(—ut) sin + 

where Rq is the initial radius in polar coordinates. Then 

r{t) = Ro\J exp(2ut) cos^ + "j) + exp(—2ut) sin^ ^ . 


AJ\ ’ 


(81) 


Differentiating this with respect to time and using equation ([7^ shows that this is a solution to (1771) . 

One can of course calculate transition probabilities for this model as in subsection l4.3l In particular, 
by equation ([551) . the probability that the oscillator stays in the ground state after time t is 


\{h^\Ut\ho)\^ = 


dpdg 

h 


dp'dq' {\ho){ho\) (p, q)Pw (p, q, t\p',q',0){\ho) (/iq I) {?', q ') 


(82) 


where 


dpdq = hdx dy = hR dR dcp , 

and, recalling equation (jS]), 

{\ho){ho\){p,q) = 2exp (-(x^ + y^)) = 2exp{-R^). 

Now 

Pw{p,q,t\p',q',0) = 6{p - p{t))5{q - q{t)) = -^hiR - r{t)6{(p - 9{t)). 

where, {r{t),9(t)) are the solutions for the parametric oscillator in polar coordinates with initial con¬ 
ditions {R',4>'). The Weyl transform of \ho){ho\ depends only on the radial coordinate and for this 
motion, using (ISTI) . 


r^(t) = R'^ exp(2ut) cos^ + T) + exp(—2ut) sin^ + T 


vr 


4/1 
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With this information (j82p simplifies to 


\{ho\Ut\ho)\^ = 

r d(/)' 
J-n TT 


exp 


—( exp(2Mt) cos^ ( + T ) + exp(—2rif) sin^ ( (f)' + 


TT 


,;„2 / 1 / 


TT 


1 + exp(2ttf) cos^ (0' + ^) + exp(—2rif) sin^ ((/>' + j) 
dcj)' 1 


/_^ vr sinh(2ttf) (cos(2(/)) + coth(ttt)) 


To evaluate the last integral here we can, as before in subsection l5.2l let z = exp(2i(;/!)) and integrate with 
respect to z twice around the unit circle in the complex z-plane, so giving (with some manipulation) 


\{ho\Ut\ho)\^ 


dz 


sinh(2nf) Jc 27ri (z — z+){z — z“) ’ 


where z^ = — tanh(nf/2) and z = — coth(nf/2). Examination shows that z^ lies within the unit 
circle and 2 ;“ without, so that, by standard means (and after a little simplification) 


\{homho)\^ 


1 

cosh(uf) 


(83) 


In this case, then, the probability that the oscillator stays in the ground state decays exponentially at 
long times. 


7 Discussion 

Ultimately, quantum mechanics calculates expectations and probabilities through the density matrix 
p{t) which always involves products of wave functions. The Wigner function pu}{p,q;t) was originally 
devised m as a scheme for using p to extract from the density matrix quantum corrections to clas¬ 
sical physics. In precise terms, the Wigner function is the Weyl transform of p/h, that is to say 
Pw < —^ p{t)/h. For systems govererned by generally time-dependent quadratic Hamiltonians the time- 
development of Pw{p,q]t) is, at its core, purely classical, with the fully quantum result to be be found 
through equations (fTSl) and (f2T)l . 

Here, we have drawn upon this quantum formulation to calculate exactly several expectations and 
probabilities for the harmonic oscillator when acted upon by an external white noise force. To do this 
one must calculate the ensemble average of p and not of the wave functions themselves. As one might 
expect, when the oscillator is initially in its ground state, over time the random force pumps it up to 
higher energy levels, statistically speaking, according to equation (l50l) . 

The ultimate randomizing effect of the stochastic force on the Weyl quantization of the oscillator’s 
angle cj) <—> (p in the phase plane, and indeed of all functions $ <—> $(</>), is reasonable: over time, 
whatever its initial state /5(0), the distribution of angle is ‘smeared’ to uniformity over the interval 
\—tt,tt). By the same token one might expect that the same force would, in time, pump up the 
oscillator’s energy. The quantum statement of this is given by ()63p . which shows, for any initial state, 
that large times 

Tr(p(f)exp(-/3.H)) ~ ^ , 

or, expanding equation ([63|) in powers of /? gives to first order, that 

Tr(/5(t)i?) —)• ^Nujt 
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where the strength of the random force is characterized by the dimensionless parameter y/N, equation 

dan). 

We also considered, in section [H an oscillator subjected not to any external force, but to a small 
time-dependent perturbation in its fundamental frequency. When the dependence is slow compared to 
the oscillator’s natural frequency—the adiabatic limit—the resultant effect on probabilities is small. 
But when, as for a parametric oscillator, this sinusoidal time-dependence has about twice the oscillator’s 
natural frequency, though the perturbation may be weak, the effects are profound. Equation (l83l) shows 
for instance that the probability that the oscillator stays in its ground state decays, for long times, 
exponentially. 
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8 Figures 



Figure 1: The expectation of phase (equation (|58l) with ^{(p) = 4>) with = 1, and cot ranging from 
0 to 20. 
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